Non-equilibrium Effects in the Thermal Switching of Underdamped Josephson Junctions 
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We study the thermal escape problem in the low damping limit. We find that finiteness of the barrier is crucial 
for explaining the thermal activation results. In this regime low barrier non-equilibrium corrections to the usual 
theories become necessary. We propose a simple theoretical extension accounting for these non-equilibrium 
processes which agrees numerical results. We apply our theory to the understanding of switching current curves 
in underdamped Josephson junctions. 
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In 1940 Kramers derived his famous formulas describing 
rates in chemical reactions 111]. The theoretical framework 
for his calculation was the escape of a Brownian particle over 
a potential barrier Far from being a particular case, noise- 
activated escape is applicable in a wide number of problems 
in science, going from biology to quantum information pro- 
cessing 1:2]. Due to the many fields involved, intense activ- 
ity emerged in the subject proposing better theories for this, 
nowadays, old problem 12|-|4[|. 

In particular, studies on thermal switching in Josephson 
junctions (JJ) benefits from this effort H- 12|. Experimental 
results are affected by thermal fluctuations and measurements 
in the lab allow to predict junctions parameters by fitting the 
switching with available expressions. Also, some fundamen- 
tal issues as the quantum-classical transition have been ad- 
dressed by means of rates measurements H-ll, 13]. It is 
clear, that such a measurements need to be compared with ap- 
propriate theoretical results. Needless to say the exact formula 
does not exist and many theories are available in the literature, 
which starting from the Kramers seminal work, cover differ- 
ent set of parameters ^ 0] ■ 

In a recent numerical work lfl4l], for very low values of the 
damping parameter, it has been found a significant deviation 
of the JJ switching current from the expected result. Here we 
present a theory that is able to give account for the observed 
deviation. Moreover, we predict that this effect will appear in 
any biased system where the damping over force ramp ratio is 
not large. In such a case, the usual theories are not suitable, 
and, as we show, it is needed to include non-equilibrium ef- 
fects and finite barrier correction in a full description of the 
problem. 

To be definite, the dynamics for the phase difference in the 
junction is usually described by the, so called, resistively and 
capacitively shunted junction (RCSJ) model, which is equiv- 
alent to the more general problem of a Brownian particle in a 
metastable potential: 



mx + mjx - 



dV . , , 



(1) 



consider white thermal noise, (^(f)} = and (^(f)^(f')) = 
2mykBTh{t~t'). 

For moderate to low values of the damping parameter there 
exists a temperature dependent critical current (force) for 
the system to switch from a superconducting or locked state 
((i)=0) to a resistive or running one {my{x) = /). Such a situ- 
ation corresponds to the problem of escape from a metastable 
well. In switching current experiments many current-voltage 
(force-velocity) curves are performed to obtain the switching 
current probability distribution function, P{I). The measured 
P{I) can be directly related to the thermal activation rate ^ 
and experimental results can be understood in terms of such 
parameter. 

For very weak damping, the Kramers result for the activa- 
tion rate is r^^^ = (Y7f,/ytBr)(co„/27t)e"'^^/^«^. There we rec- 
ognize the transition-state-theory result multiplied by a pref- 
actor valid in the very low damping regime. For our system 
the action at the barrier Jb is usually approached by the cubic 
potential result Jy = 7.2AU / (On- Then: 



(2) 



where the potential V{x) = Vo{l — cosx) — Ix and ^(f) is the 
stochastic force describing the thermal fluctuations. Here we 



This equation shows that the rate scales linearly with the 
damping and depends only on the barrier height over ther- 
mal energy factor, AU /ksT . This expression is only valid in 
the low damping and infinite barrier limit {yJb/k-aT <^ 1 and 
MJ/k^T > 1). 

Many theories have extended the Kramers result to the 
moderate-to-small damping regime following the 

infinite barrier approximation. Given its simplicity, the re- 
sult of Biittiker et al. (BHL) [fisll has been usually applied 
in the JJ Uterature. Finite barrier corrections have been stud- 
ied in ]19, 2^]. More recently, Drozdov and Hayashi (DH) 
proposed a new theory which is not perturbative in the barrier 
height i2lil . 

We are interested in the dynamics of the system in the low 
damping limit. In this limit the coupling to the bath is very 
weak and the time to reach thermal equilibrium very long 
1/y). This fact has important consequences: For biased 
systems, escape occurs at very low values of the At/ /ksT ratio 
and junctions may escape before thermal equilibrium is estab- 
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FIG. 1 : (color online) Escape rate divided by the Kramers low damp- 
ing result, Eq. [21 as a function of MJ /k^T. Figure (a): Dots are 
numerical results and the corresponding color solid lines are the the- 
oretical prediction by DH theory 12211 . The black line stands for the 
vanishing damping rate formula, Eq. Figures (b) and (c) com- 
pare different theoretical results. 



lished and thus non-equilibrium effects dominate the process. 
In order to study such effects we need first to know the im- 
portance of finite barrier effects in particle activation problem 
at low damping and take into account the average energy of 
junctions before each switching event. 

Escape at small barrier — We will show here that small 
barrier effects are very important in the low damping case, the 
convergence to the infinite barrier result is very slow and the 
DH theory reproduces the numerical results at any barrier 

We have numerically integrated the Langevin equation ([Hi 
of the system for different values of damping and barrier 
height. In our simulations we have computed the mean time 
for the system to first reach the potential barrier For low val- 
ues of the damping such mean time (the first-passage-time 
problem) corresponds to the inverse of the escape rate. Ac- 
cording to theory, simulations are started with particles placed 
in the metastable potential well and zero velocity. Some issues 
about the initial conditions problem will be addressed below. 
At any point the numerical result is obtained from 10^ escap- 
ing events. We show results for Vo ~ 0.31, m ~ 0.35 and dif- 
ferent values of F, damping and temperature. The results are 
summarized in Figs.[T]and[2|where we plot the activation rate 
as a function of barrier and damping respectively and compare 
to some existing theories. 

Figure [T] shows the rate dependence on the barrier for dif- 



FIG. 2: (color online) Damping dependence for the escape rates. 
The points are numerical results and the corresponding solid lines 
account for DH theory 1I22I1 . The horizontal black line is the vanish- 
ing damping limit Eq. ([Sj. Top: curves for At/ /k^T = 3. Bottom: 
AU /ksT = 10. The inset compares DH to BHL (green lines) results. 



ferent values of damping and temperature. In order to see de- 
viations from the Kramers low damping result we divide the 
obtained rates by Eq. We recall that ^ is obtained as- 
suming weak damping and high barrier For comparison, we 
also plot the exact result for arbitrary barrier in the limit of 
vanishing damping jstl, r{j^O) = rnrB with 



-mj) 



J 2n 



-1 



(3) 

where Jj, is the action at the barrier and p = l/ksT . 

Remarkably, the approach between both results is slow, 
meaning that the high barrier approximation is accurate only 
at very high barrier values indeed [Fig. [It a)]- For instance 
rHTB/rKLD 0.72 for MJ /keT = 5 and 0.85 for AU/kBT = 
10. As a consequence, all theories which try to extend 
the Kramers result to the moderate-to-small damping re- 
gion is, 15-18] also fail at low damping values unless very 
large barriers are considered. 

We know about two main attempts to include finite bar- 
rier effects in this limit: the first one is due to Melnikov 
(MFB) iH] and fails at small barriers, see Fig. [lib). The 
second one was proposed by Drozdov and Hayashi (DH) for 
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FIG. 3: (color online) Activation rate vs barrier for three initial condi- 
tions: Particles are at the bottom of the metastable well with zero ve- 
locity (blue points), v = ^Jk^T /m (green points) or v = —^/k^Tjm 
(brown points). Lines stand for theoretical predictions of Eq. lO. 



moderate-to-small damping and arbitrary barrier H21L 12211 . As 
it can be seen in Figs. [T]and|2] the DH theory recovers the 
vanishing damping limit and explain our numerical results in 
the whole small damping region. In Figslljb) and (c) we show 
BHL result for a = 1 . Other extensions of the Kramers low 
damping result to the moderate-to-low damping regime give a 
result quite similar to BHL. 

In Fig. 12] we plot the rate damping dependence for two dif- 
ferent barriers. Apart from the agreement to the DH theory, 
we check the convergence to the very weak limit given by ^ 
[Cf. horizontal line in both figures]. In the inset we also show 
results for the BHL. We see that the finite barrier corrections 
become less important by increasing the damping. 

To complete our discussion we notice that by decreasing 
the damping the curves at different temperatures become the 
same, that is the rate depends only on the ratio AU / k^T . This 
can be understood by noticing that most of the contribution 
in the integrals in Eq. ([3j comes from the bottom of the po- 
tential. If the action inside the well is approximated by the 
corresponding one for an harmonic potential with the same 
bottom frequency, E = {(£io/2n)J, with coo — 9^V(x„„„), the 
rate Q can be written as. 



r-HA 



dx e^'^ / dy — 

Jx y 



(4) 



which clearly depends only on the ratio AU /k^T . Let us em- 
phasize that, besides its simplicity, the above equation is an 
excellent approximation to (O. In fact, plotting both (|4|i and 
(|3]l in Fig.[T] they cannot be distinguished one from the other. 

Initial Conditions dependence — We address now the influ- 
ence of the initial condition for the energy on the escape rate 
results. Thermal escape at low damping is an energy diffusion 
problem. Escape occurs as soon as thermal fluctuations pro- 
vide a particle energy enough to overcome the barrier. This 



time depends on the value of the particle initial energy. 

Up to know, to compare our simulations to theory we as- 
sumed that the particle starts at the bottom of the metastable 
well with zero velocity. From the experimental point of view 
this assumption may fail. Thermal fluctuations not only pro- 
vide energy enough to surmount the barrier but also kinetic 
energy at the bottom. In order to study the importance of this 
issue, in Fig. |3]we plot the rates with two initial conditions, 
V = ±^ ksT / m and compare to the one with zero velocity. 
As expected, we see that for small barriers initial kinetic en- 
ergy speed up the escape times. 

When particles are placed with zero velocity at the bottom 
of the well, the activation time is r^^ . However, if particles 
have extra initial energy £,„ the escape time is given by ' — x 
where x is the activation time up to this extra energy, which 
can be computed at low damping from Eq.© replacing AU 
by £■,„. Putting all together we generalize the rate formulas as. 



1 



(5) 



This equation shows that as soon as t(£',„) ^ r^' the ini- 
tial conditions problem affect the escape rates. In Fig. |3] 
where £„, = k^T /2, this correction becomes important for 
AU /k^T < 2. If AU < Ei„ the passage time is almost a de- 
terministic process which depends on the initial position and 
velocity of the system. Figure[3]illustrates this effect and con- 
firms our theoretical prediction. 

Switching current — In a typical JJ experiment the probabil- 
ity distribution function of the junction switching current P{I) 
is measured performing many current-voltage curves where 
current is continuously increased at a given rate. From these 
results the mean switching current /sw and its standard devia- 
tion can be trivially computed. Such P{I) can be easily related 
to the escape rates r{I) as jsl] 



-1 



P{u)du 



(6) 



Alike, escape rates can be computed from measured P{I). 

Figure Ufa) shows our numerical results for the average 
switching current and compares them to theoretical predic- 
tions. We integrate Eq. ([T]) for an ensemble of thermalized 
junctions. Current is increased at a given ramp and switch- 
ing events are recorded. As expected BHL based predictions 
fail in the very low damping regime. However, surprisingly, 
also DH is unable to explain our numerical results, which lie 
in between both theories. This is due to the competition be- 
tween the equilibrium time of the system, given by y^', and 
the time order for the change of the current, given by the in- 
verse of the current ramp. Thus, switching in the very low 
damping regime is a non-equilibrium process. The coupling 
to the external bath is so weak that other junctions are not able 
to reach the thermal energy before switching. Thus, junctions 
escape in an evaporative cooling way where more energetic 
junctions switch first and the ensemble is effectively cooled. 
This picture is confirmed in Fig.|4|c) where for a given damp- 
ing we show the mean energy for the trapped junctions as a 
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FIG. 4: (color online) (a) Average switching current at different val- 
ues of the damping for T = 0.01 and current ramp / = 3.33 x 10^^ 
(red) and / = 10^^ (blue). Solid symbols are for numerical sim- 
ulations, dashed lines for predictions based on BHL theory, dot- 
ted lines for DH theory and solid lines for our theory, Eq. ^ with 
Ein = kgT /2. We also plot (right axis) the value of the barrier at the 
mean switching current, (b) Normalized current versus normalized 
damping over ramp ratio (23ll . (c) Mean energy divided by kgT of 
particles in the well (brown) and fraction of particles which have es- 
caped (green) as a function of the applied current {T = 0.01, Y= 10^^ 
and / = 3.33 X 10"^). 



function of the cuiTent and the fraction of particles which have 
switched. 

We also see in Fig. Uc) that particles escape with an ini- 
tial energy which goes from E\a ~ KqT Xo zero when current 
is increased. The simplest way to introduce this fact in the 
theory is to assume an average value of E\n = KqT /2 and use 
our Eq. (|5]). Figure HJa) shows that in this way we are able 
to reproduce quite accurately the numerical results. This cor- 
rection turns out to be important when the average barrier at 
the switching cuiTent is of the order of the thermal energy. 
See that in the figure it is also plotted the value of the barrier 
at the mean switching current (open symbols). Finally, using 
Eqs. (HI, ^ and ^ it can be seen that in this region of the 
parameter space, the results depend on the y// ratio, as con- 
firmed in Fig. lUb). Therefore our theory allows to estimate 
the values for y// where non equilibrium corrections are nec- 
essary. 

Although presented in the framework of the JJ switching 
current measurements, our results are further more general 
and apply to any experiment where an activation rate is mea- 
sured as a function of an external parameter which can be con- 
trolled at will. An important issue to study, it is the influence 



of the observed competition between two different time scales 
on results for biased systems at higher values of the damp- 
ing and thus transfer our theoretical scheme from the energy- 
diffusion regime to the phase-diffusion one. This is the typ- 
ical case for many of the current biological-physics experi- 
ments ll24il25ll 
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